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The liquid-gas phase transition in hot neutron-rich nuclear matter is investigated within a self- 
consistent thermal model using an isospin and momentum dependent interaction (MDI) constrained 
by the isospin diffusion data in heavy-ion collisions, a momentum- independent interaction (MID), 
and an isoscalar momentum-dependent interaction (eMDYI). The boundary of the phase-coexistence 
region is shown to be sensitive to the density dependence of the nuclear symmetry energy with a 
softer symmetry energy giving a higher critical pressure and a larger area of phase-coexistence region. 
Compared with the momentum-independent MID interaction, the isospin and momentum-dependent 
MDI interaction is found to increase the critical pressure and enlarge the area of phase-coexistence 
region. For the isoscalar momentum-dependent eMDYI interaction, a limiting pressure above which 
the liquid-gas phase transition cannot take place has been found and it is shown to be sensitive to 
the stiffness of the symmetry energy. 
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The van der Waals behavior of the nucleon-nucleon in- 
teraction is expected to lead to the so-called liquid-gas 
(LG) phase transition in nuclear matter. Since the early 
work, see, e.g., Refs. [l|, H, H, many investigations 
have been carried out to explore properties of the nu- 
clear LG phase transition both experimentally and the- 
oretically over the last three decades. For a recent re- 
view, see, e.g., Refs. d, d, 0- Most of these studies 
focused on investigating features of the LG phase tran- 
sition in symmetric nuclear matter. Recent progress in 
experiments with radioactive beams provided us a great 
opportunity to explore more extensively the LG phase 
transition in isospin-asymmetric nuclear matter. New 
features of the LG phase transition in asymmetric nu- 
clear matter are expected. In particular, in an asym- 
metric nuclear matter, there are two components of pro- 
tons and neutrons, and two conserved charges of baryon 
number and the third component of isospin, and the 
LG phase transition was suggested to be of second or- 
der [8|. This suggestion together with the need to un- 
derstand better properties of asymmetric nuclear mat- 
ter relevant for both nuclear physics and astrophysics 
have stimulated a lot of work recently, see, e.g., Refs. 
S E, [H m, [11 m, [S [M [13, WhUe signifi- 

cant progress has been made recently, many interesting 
questions about the nature of the LG phase transition in 
asymmetric nuclear matter remain open. Some of these 
questions can be traced back to our poor understanding 
about the isovector nuclear interaction and the de nsity 
dependence of the nuclear symmetry energy 0, [H, [20| . 
Fortunately, recent analyses of the isospin diffusion data 
in heavy-ion reactions have allowed us to put a stringent 
constraint on the symmetry energy of neutron-rich mat- 
ter at sub- normal densities [ll], [52, It is therefore 
interesting to investigate how the constrained symme- 
try energy may allow us to better understand the LG 
phase transition in asymmetric nuclear matter. More- 



over, both the isovector (i.e., the nuclear symmetry po- 
tential) and isoscalar parts of the single nucleon potential 
should be momentum dependent due to the non-locality 
of nucleon-nucleon interaction and the Pauli exchange 
effects in many-fermion systems. However, effects of 
the momentum-dependent interactions on the LG phase 
transition in asymmetric nuclear matter have received so 
far little theoretical attention. A timely remedy of the 
situation is imperative given the experimental studies un- 
derway or planned at various radioactive beam facilities. 

In the present work, we study effects of isospin and mo- 
mentum dependent interactions on the LG phase tran- 
sition in hot neutron-rich nuclear matter within a self- 
consistent thermal model using three different interac- 
tions. The first one is the isospin and momentum de- 
pendent MDI interaction constrained by the isospin dif- 
fusion data in heavy-ion collisions. The second one is a 
momentum-independent interaction (MID) which leads 
to a fully momentum independent single nucleon po- 
tential, and the third one is an isoscalar momentum- 
dependent interaction (eMDYI) in which the isoscalar 
part of the single nucleon potential is momentum depen- 
dent but the isovector part of the single nucleon potential 
is momentum independent. 

In the isospin and momentum-dependent MDI inter- 
action, the potential energy density Vmdi(p, 2^,(5) of a 
thermally equilibrated asymmetric nuclear matter at to- 
tal density p, temperature T and isospin asymmetry d is 
expressed as follows [H, [23| , 
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In the mean field approximation, Eq. ((T]) leads to the 
following single particle potential for a nucleon with mo- 
mentum p and isospin r in the thermally equilibrated 
asymmetric nuclear matter (22l . [23 | 
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In the above r = 1/2 (—1/2) for neutrons (protons) and 
fri^jP) is the phase space distribution function at coor- 
dinate r and momentum p. The detailed values of the 
parameters g, A„ ( x), Ai (x), B, Cr.r, Cr.-r and A can be 
found in Refs. [22,|2J] and have been assumed to be tem- 
perature independent here. The isospin and momentum- 
dependent MDI interaction gives the binding energy per 
nucleon of —16 MeV, incompressibility Kg = 211 MeV 
and the symmetry energy of 31.6 MeV for cold symmet- 
ric nuclear matter at saturation density pa = 0.16 fm^"^ 
[23 |. The different x values in the MDI interaction are 
introduced to vary the density dependence of the nuclear 
symmetry energy while keeping other properties of the 
nuclear equation of state fixed [l^l- We note that the 
MDI interaction has been extensively used in the trans- 
port model for studying isospin effects in intermediate 
ener gy heavy-ion collisions induced by neutron-rich nu- 
clei [ii, H, m, m, [13, m, [H [131. in particular, the 
isospin diffusion data from NSCL/MSU have constrained 
the value of x to be between and —1 for nuclear mat- 
ter densities less than about 1.2po [13, [13, we will thus 
in the present work consider the two values of a; = 
and X — —1. We note that the zero-temperature sym- 
metry energy for the MDI interaction with x — and 
— 1 can be parameterized, respectively, as 31.6(p//9o)*''®^ 
MeV and 31.6(p/po)^ °^ MeV [H, and thus x = gives 
a softer symmetry energy while x = —1 gives a stiffer 
symmetry energy. 

In the momentum-independent MID interaction, the 
potential energy density Vuiuip, <5) of a thermally equili- 
brated asymmetric nuclear matter at total density p and 
isospin asymmetry S can be written as 



Vuiu {p,S)^-^ + -r^^+ pEi;l {p, x)6' . (3) 
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The parameters a, (3 and 7 are determined by the in- 
compressibility Kq of cold symmetric nuclear matter at 
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and Kq is again set to be 211 MeV as in the MDI inter- 
action. To fit the MDI interaction at zero temperature, 
the potential part of the symmetry energy i?f (p, x) is 
parameterized by (2^ 

Erylip, x) = F{x)^ + [18.6 - Fix)] (-^)G(-) (7) 
Po Po 

with F{x = 0) 129.981 MeV, G{x = 0) = 1.059, 
F{x = -1) 3.673 MeV, and G{x = -1) = 1.569. 
We note that the MID interaction reproduces very well 
the EOS of isospin-asymmetric nuclear matter with the 
MDI interaction at zero temperature for both x = and 
X — —1. The single nucleon potential in the MID inter- 
action can be directly obtained as 
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Therefore, the single nucleon potential in the MID inter- 
action is fully momentum-independent. It also leads to 
the fact that the potential energy density and the single 
nucleon potential in the MID interaction are independent 
of the temperature. 

The momentum-dependent part in the MDI interac- 
tion is also isospin dependent while the MID interac- 
tion is fully momentum independent. In order to see 
the effect of the momentum dependence of the isovec- 
tor part of the single nucleon potential (nuclear symme- 
try potential), we can introduce an isoscalar momentum- 
dependent interaction, called extended MDYI (eMDYI) 
interaction since it has the same functional form as the 
well-known MDYI interaction for symmetric nuclear mat- 
ter [m. In the eMDYI interaction, the potential energy 
density I^mdyKp, T, 5) of a thermally equilibrated asym- 
metric nuclear matter at total density p, temperature T 
and isospin asymmetry 5 can be written as 
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Here fo{r,p) is the phase space distribution function of 
symmetric nuclear matter at total density p and temper- 
ature T. EPy!^{p,x) has the same expression as Eq. 
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have the same values as in the MDI interaction, so that 

the eMDYI interaction gives the same EOS of asymmet- 
ric nuclear matter as the MDI interaction at zero tem- 
perature for both X = and x = —1. The single nucleon 
potential in the eMDYI interaction can be obtained as 

U,MUYi{p,T,S,p,T)^U''{p,T,p) + U'''y{p,6,T), (11) 
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and If^^^lp, S, t) is the same as Eq. ^ which implies that 
the symmetry potential is identical for the eMDYI and 
MID interactions. Therefore, in the eMDYI interaction, 
the isoscalar part of the single nucleon potential is mo- 
mentum dependent but the nuclear symmetry potential 
is not. For symmetric nuclear matter, the single nucleon 
potential in the eMDYI interaction is exactly the same 
as that in the MDI interaction. We note that the same 
strategy has been used to study the momentum depen- 
dence effects in heavy-ion collisions in a previous work 

At zero temperature, fT{r,p) = ■p-6(p/(r) — p) and all 
the integrals in above expressions can be calculated ana- 
lytically, while at a finite temperature T, the phase space 
distribution function becomes the Fermi distribution 
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where pr is the proton or neutron chemical potential and 
can be determined from 



(14) 



In the above, is the proton or neutron mass and Ut is 
the proton or neutron single nucleon potential in different 
interactions. From a self-consistency iteration scheme 
[31I [3^ , the chemical potential pr and the distribution 
function fri^jp) can be determined numerically. 

From the chemical potential Pt- and the distribution 
function fr{r,p)^ the energy per nucleon E{p,T,S) can 
be obtained as 



E{p,T,S) = - 
P 
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Furthermore, we can obtain the entropy per nucleon 
St [p, T, 5) as 
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with the occupation probability 
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Finally, the pressure P{p, T, S) can be calculated from 
the thermodynamic relation 
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With the above theoretical models, we can now study 
the LG phase transition in hot asymmetric nuclear mat- 
ter. The phase coexistence is governed by the Gibbs con- 
ditions and for the asymmetric nuclear matter two-phase 
coexistence equations are 

pf{T, pf) = pfiT, pf), {i^n and p) (19) 
Pt{T, ) = Pf^{T, pf), {i = nov p) (20) 

where L and G stand for liquid phase and gas phase, 
respectively. The chemical stability condition is given by 
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The Gibbs conditions (fT9| and ([20|) for phase equilibrium 
require equal pressures and chemical potentials for two 
phases with different concentrations. For a fixed pres- 
sure, the two solutions thus form the edges of a rectangle 
in the proton and neutron chemical potential isobars as 
a function of isospin asymmetry S and can be found by 
means of the geometrical construction method [1, [ll| . 

The solid curves shown in Fig. [1] are the proton and 
neutron chemical potential isobars as a function of the 
isospin asymmetry 6 at a fixed temperature T — 10 MeV 
and pressure P = 0.090 MeV/fm'^ by using the MDI and 
MID interactions with x = and x = —1. The result- 
ing rectangles from the geometrical construction are also 
shown by dotted lines in Fig. [1] For each interaction, 
the two different values of 6 correspond to two differ- 
ent phases with different densities and the lower den- 
sity phase (with larger S value) defines a gas phase while 
the higher density phase (with smaller 6 value) defines 
a liquid phase. Collecting all such pairs of S{T, P) and 
(5'(T, P) thus forms the binodal surface. From Fig.[Tl one 
can see that different interactions give different shapes 
for the chemical potential isobar. When the pressure in- 
creases and approaches the critical pressure Pc , an inflec- 
tion point will appear for proton and neutron chemical 
potential isobars, i.e., 
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Above the critical pressure, the chemical potential of neu- 
trons (protons) increases (decreases) monotonically with 



4 









1 1 1 1 1 1 1 




































MDI with X = 




MIDwithx = 






P= 


0.090 MeV/fm' 




P=0.090 MeV/fm' 


\ \ 


p= 


0.265 MeV/fm' 




" P=0.230 MeV/fm' 




— 1 {- 












. ^^^^^^^^...^^^ 




MDI with X = -1 




'-^~--<^,^ 
MID with X = -1 






P= 


0.090 MeV/fm' 




" P=0.090 MeV/fm 






p= 

. 


0.195 MeV/fm' 




" P=0.154MeV/fm 









0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 
8 

FIG. 1: (Color online) The chemical potential isobar as a 
function of the isospin asymmetry 5 at T = 10 MeV in the 
MDI and MID interactions with x — and x — —1. The geo- 
metrical construction used to obtain the isospin asymmetries 
and chemical potentials in the two coexisting phases is also 
shown. 



d and the chemical instability disappears. In Fig. [l] we 
also show the chemical potential isobar at the critical 
pressure by the dashed curves. At the critical pressure, 
the rectangle is degenerated to a line vertical to the 6 
axis as shown by dash-dotted lines in Fig. [TJ The values 
of the critical pressure are 0.265, 0.230, 0.195 and 0.154 
MeV/fm^ for the MDI interaction with a; = 0, MID in- 
teraction with a; = 0, MDI interaction with x = —I and 
MID interaction with x — —1, respectively. It is inter- 
esting to see that the different interactions give different 
values of the critical pressure. Especially, the stiffer sym- 
metry energy {x = —1) gives a smaller critical pressure 
while the inclusion of the momentum dependence in the 
interaction (MDI) gives a larger critical pressure. 

In Fig. [21 we show the section of the binodal surface 
at T = 10 MeV for the MDI and MID interactions with 
a; = and x — —1. On the left side of the binodal 
surface there only exists a liquid phase and on the right 
side only a gas phase exists. In the region of "filet 
mignon" is the coexistence phase of liquid phase and 
gas phase. Interestingly, we can see from Fig. [2] that the 
stiffer symmetry energy (x = —1) significantly lowers the 
critical point (CP) and makes the maximal asymmetry 
(MA) a little smaller. Meanwhile, the momentum 
dependence in the interaction (MDI) lifts the CP in a 
larger amount, while it seems to have no effects on the 
MA point. In addition, just as expected, the value of x 
does not affect the equal concentration (EC) point while 
the momentum dependence lifts it slightly (by about 
0.005 MeV/fm'^). These features clearly indicate that 
the critical pressure and the area of phase-coexistence 
region in hot asymmetric nuclear matter is very sensitive 
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FIG. 2: (Color online) The section of binodal surface at T = 
10 MeV in the MDI and MID interactions with x = and x = 
— 1. The critical point (CP), the points of equal concentration 
(EC) and maximal asymmetry (MA) are also indicated. 

to the stiffness of the symmetry energy with a softer 
symmetry energy giving a higher critical pressure and a 
larger area of phase-coexistence region. Meanwhile, the 
momentum dependence in the interaction significantly 
increases the critical pressure and enlarges the area of 
phase-coexistence region. 




FIG. 3: (Color online) The chemical potential isobar as a 
function of isospin asymmetry 5 at T = 10 MeV in the eMDYI 
interaction with a:: = and x = —1. 



We now turn to the case of the eMDYI interaction. In 
the eMDYI interaction, the resulting single nucleon po- 
tential is momentum dependent but its momentum de- 
pendence is isospin-independent. Comparing the results 
with those of the MDI interaction, we can extract infor- 
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mation about the effects of the momentum dependence of 
the symmetry potential while the effects of the momen- 
tum dependence of the isoscalar part of the single nu- 
cleon potential can be studied by comparing the results 
with those of the MID interaction. Shown in Fig. [3] is 
the chemical potential isobar as a function of the isospin 
asymmetry (5 at T = 10 MeV by using the eMDYI in- 
teraction with a; = and x — —1. Compared with the 
results from the MDI and MID interactions, the main dif- 
ference is that the left (and right) extrema of /i„ and fip 
do not correspond to a same 6 but they do for the MDI 
and MID interactions as shown in Fig. [T] In particular, 
for the eMDYI interaction, the chemical potentials of pro- 
tons and neutrons are seen to exhibit asynchronous varia- 
tion with pressure, i.e., the chemical potential of neutrons 
increase more rapidly with pressure than that of protons. 
This asynchronous variation is uniquely determined by 
the special momentum dependence in the eMDYI inter- 
action within the present self-consistent thermal model. 
Actually, it is this asynchronous variation that leads to 
the fact that the left (and right) extrema of /i„ and fip 
correspond to different values of S. 

At lower pressures, for example, P — 0.090 MeV/fm"^ 
as shown in Fig. [3] (a), the rectangle can be accurately 
constructed and thus the Gibbs conditions and ([^0]) 
have two solutions. Due to the asynchronous variation 
of Hn and /Zp with pressure, we will get a limiting pres- 
sure Piim above which no rectangle can be constructed 
and the coexistence equations (|19p and pp]) have no so- 
lution. Fig. [3] (b) shows the case at the limiting pressure 
with Pii„i = 0.205 and 0.175 MeV/fm^ for a; = and 
X = —I, respectively. In this limit case, we note that for 
a; = the left edge of the rectangle actually corresponds 
to the left extremum of /ip and the pair (5i.o — 0.600 and 
'52.0 = 0.750 form the two edges of the rectangle (the end 
of the binodal surface) while for x = — 1 it corresponds 
to the left extremum of fin and the pair (5i._i — 0.480 
and 62,-1 = 0.811 form the two edges of the rectan- 
gle. The ^3,0 = 0.720 and S3-1 = 0.700 indicated in 
Fig. [3] (b) correspond to the maximum of /ip for a; = 
and X — —1, respectively. With increasing pressure, 
namely, at P = 0.260 and 0.230 MeV/fm^ for a; = 
and X = — 1, respectively, fi„ passes through an inflec- 
tion point while /ip still has a chemically unstable region 
and this case is shown in Fig. [3] (c). When the pressure 
is further increased to P = 0.305 and 0.250 MeV/fm^ for 
a; = and a; = — 1, respectively, as shown in Fig. [3] (d), 
Hp passes through an inflection point while /x„ increases 
monotonically with 6. These features indicate that the 
asynchronous variation of fin and fip with pressure also 
depends on the value of x. 

Fig. d] displays the section of the binodal surface at 
T = 10 MeV by using the eMDYI interaction with a; = 
and a; = —1. We can see that the curve is cut off 
at the limiting pressure with Piim = 0.205 and 0.175 
MeV/fm^ for a; = and x = —1, respectively. It is 
interesting to see that for a; = — 1 there is an bend- 
ing point at Srnax.-i — 0.823, but for x = there is 
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FIG. 4: (Color online) The section of binodal surface at T = 
10 MeV in the eMDYI interaction with a:: = and x = — 1. 
Piini represents the limiting pressure and 5ij 's are discussed 
in the text. 



no such point. For a; = — 1 the situation is similar to 
that found in Ref. till] using different models and the 
isospin asymmetry 8 can be divided into four regions: 
[0, -1,(^3,-1], [(^3, -1,(52,-1] and [(52,-1,(5 

In the first region [0, (5i._i], the system begins at gas 
phase, experiences a LG phase transition and ends at 
liquid phase. In the second region [(5i._i, (53._i], it begins 
at gas phase, enters a two-phase region and becomes un- 
stable at the limiting pressure, as chemical instability 
condition Eq. pT|) is destroyed in this region. In the 
third region [(53, -1, (52, -1], the system will end in a stable 
phase as Eq. (pij) is satisfied in this region. In the fourth 
region [(^2,-1, (5maa;,-i], the system enters and leaves the 
two-phase region on the same branch, so it remains in 
the gas phase at the end. For a; = 0, the situation is 
simpler and b can be divided into three regions: [0, (5i,o], 
[<5i,o, <53_o] and [(53^_i, (52,o]- The discussion is similar to 
that of the first, second and third regions in the case of 
X = -1. 

From Fig. 31 we can also see that the limiting pressure 
and the area of phase-coexistence region are still sensitive 
to the stiffness of the symmetry energy with a softer sym- 
metry energy (a: — 0) giving a higher limit pressure and 
a larger area of phase-coexistence region. Comparing the 
results of the MDI and MID interactions shown in Fig.[2l 
we can see that for pressures lower than the limiting pres- 
sure, the binodal surface from the eMDYI interaction is 
similar to that from the MDI interaction. This feature 
implies that the momentum dependence of the symmetry 
potential has little influence on the LG phase transition 
in hot asymmetric nuclear matter while the momentum 
dependence of the isoscalar single nucleon potential sig- 
nificantly enlarges the area of phase-coexistence region 
for pressures lower than the limiting pressure. For pres- 
sures above the limiting pressure, the momentum depen- 
dence of both the isoscalar and isovector single nucleon 
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potentials becomes important. 

In summary, we have studied the hquid-gas phase 
transition in hot neutron-rich nuclear matter within 
a self-consistent thermal model using three differ- 
ent nuclear effective interactions, namely, the isospin 
and momentum dependent MDI interaction constrained 
by the isospin diffusion data in heavy-ion collisions, 
the momentum-independent MID interaction, and the 
isoscalar momentum-dependent eMDYI interaction. At 
zero temperature, the above three interactions give the 
same EOS for asymmetric nuclear matter. The MDI in- 
teraction is realistic, while the MID and eMDYI inter- 
actions are only used as references in order to explore 
effects of the isospin and momentum dependence on the 
liquid-gas phase transition. Since the symmetry energy 
corresponding to the MDI interaction has been already 
restricted within a narrow range by the isospin diffusion 
data in heavy- ion collisions, predictions using the MDI 
interaction thus allow us to constrain the liquid-gas phase 
boundary in asymmetric nuclear matter. 

Comparing calculations with the three interactions, we 
find that the boundary of the phase-coexistence region is 
very sensitive to the density dependence of the nuclear 
symmetry energy. A softer symmetry energy leads to a 
higher critical pressure and a larger area of the phase- 
coexistence region. On the other hand, compared with 
the momentum-independent MID interaction, the isospin 
and momentum-dependent MDI interaction increases the 
critical pressure and enlarge the area of phase-coexistence 



region. For the isoscalar momentum-dependent eMDYI 
interaction, there is a limiting pressure above which the 
liquid-gas phase transition cannot take place due to the 
asynchronous variations of the nucleon chemical poten- 
tials with pressure. Comparing the results of the eMDYI 
interaction with those of the MDI and MID interactions, 
we find that, for pressures lower than the limiting pres- 
sure, the momentum dependence of the symmetry poten- 
tial has little influence on the liquid-gas phase transition 
in hot asymmetric nuclear matter. While the momentum 
dependence of the isoscalar single nucleon potential sig- 
nificantly enlarges the area of phase-coexistence region. 
For pressures above the limiting pressure, the momentum 
dependence of both the isoscalar and isovector single nu- 
cleon potentials becomes important. 
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